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Vibrationally inelastic electron transport through a flexible molecular junction is in- 
vestigated. The study is based on a mechanistic model for a biphenyl molecule between 
two metal electrodes. Employing methods from electron-molecule scattering theory, which 
allow a numerically exact treatment, we study the effect of vibrational excitation on the 
transmission probability for different parameter regimes. The current-voltage characteris- 
tic is analyzed for different temperatures, based on a Landauer-type formula. Furthermore, 
the process of electron assisted tunneling between adjacent wells in the torsional poten- 
tial of the molecule is discussed and the validity of approximate methods to describe the 
transmission probability is investigated. 
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1 Introduction 

The field of molecular electronics, which may represent the ultimate limit of 
the miniaturization of electronic devices, has received much interest in the last 
decade (see for example and references therein). Although the basic idea to 

use molecules as active elements in electronic circuits is not new jS], only recently 
it became possible to study the conduction of single molecules experimentally (0J 
Him E|). In contrast to metal or semiconductor devices, in molecular electronics the 
complex structure of a single molecule may be utilized to obtain the functionality of 
the device. For example, switching behavior and negative differential conductivity 
|H1 E] have recently been demonstrated experimentally and are believed to be a 
consequence of geometry changes in the conducting molecule. 

The tremendous experimental progress has stimulated great interest in the the- 
ory and ab initio modeling of charge transport through single molecules (see for 
example [Q|2| and references therein). In particular, non-equilibrium Green's func- 
tion methods in combination with state-of-the art electronic structure calculations 
have been employed to study the conductivity of metal-molecule-metal junctions 
with fixed nuclear geometry. The dependence of the conductive properties on the 
binding geometry of the molecular bridge has also been investigated [TUinT] . 
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The role of the vibrational degrees of freedom of the molecular bridge in the con- 
duction process and and their effect on the functionality of the device are less well 
understood. The "static" influence of the internal vibrational modes has been stud- 
ied by averaging the elastic transmittance over the probability distribution of the 
vibrational degrees of freedom |12U13ll*H| . The dynamical impact of the vibrational 
degrees of freedom on the tunneling current in molecular junctions has been investi- 
gated within, e.g., nearest neighbor tight-binding models E3 El EH EH • These 
studies have demonstrated that the vibrational motion of the molecular bridge may 
result in additional (vibrational) resonance structures in the transmission probabil- 
ity which can alter the current-voltage characteristic significantly. Dynamical effects 
of nuclear motion on the conductivity have also been observed experimentally. For 
example, in experiments on electron transport through H2 molecules between two 
platinum electrodes |H] as well as Ceo molecules connected to gold electrodes [22, 
indications for an influence of the center-of-mass motion of the respective molecule 
on the conductivity have been found. Moreover, the excitation of the vibrational de- 
grees of freedom of the molecule provides a mechanism for heating of the molecular 
junction and thus is a possible source of instability |21II22| . 
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Fig. 1. Schematic illustration of the molecular junction studied. 

Recently, we have demonstrated that methods developed in the field of 
resonant electron-molecule scattering [21] can be advantageously applied to study 
vibrationally inelastic effects on electron transport. In this study, the vibrational 
degrees of freedom of the molecular bridge were described within the harmonic 
approximation. As an extension of this work, we investigate in this paper the influ- 
ence of large-amplitude motion of a flexible molecular bridge on charge transport 
through a molecular junction. Specifically, we consider a mechanistic model for 
electron conduction through a biphenyl molecule, schematically shown in Fig.E I n 
biphenyl and similar molecules (such as bypiridine jJS] or biphenyl dithiol HO]), the 
torsion of the two benzene rings is expected to have a strong effect on the conduc- 
tivity. Ratner, Datta and coworkers have considered this effect for a static nuclear 
bridge [27| as well as in the adiabatic approximation [T3|. Ab initio calculations of 
equilibrium geometries of neutral molecules and anions for several related molecules 
|S] also indicate that the torsional motion is presumably rather strongly coupled 
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to the tunneling electron. Dynamical effect of vibrations beyond the harmonic ap- 
proximation have been studied for similar systems based on molecular dynamics 
simulations [2H1, where the nuclear motion was treated classically and the elec- 
tronic structure was described by a density-functional tight-binding Hamiltonian. 
The purpose of this paper is to study the dynamical effect of the torsional motion 
on the transport through the bridge within a fully quantum mechanical treatment. 

This paper is organized as follows: In section 2.1 we introduce a generic model 
Hamiltonian to describe charge transport through a flexible molecular junction 
schematically shown in Fig. ^ The torsional angle between the two benzene rings 
is thereby treated as a dynamical variable. The model parameters are chosen in 
accordance with available ab initio data for the biphenyl molecule. The methods 
employed to calculate the inelastic transmission probability and the current for this 
model are presented in section 2.2. In section 3 we discuss the results of the model 
study for different parameter regimes and analyze the performance of approximate 
methods, in particular the purely elastic calculation of the transmission probability 
and the adiabatic-nuclei approximation. 



2.1 Model 

We describe charge transport through the molecular junction schematically de- 
picted in Fig.^in the single particle approximation as the subsequent transmission 
of single electrons from a state \<f>k), k <E L, in the conduction band of the left lead 
into state \4>k), k G R in the right lead. We assume that the conduction electron 
can cross the junction only through a 7r-orbital |</>i) associated with the benzene 
ring which is coupled directly to the left lead and through a 7r-orbital 02 ) associ- 
ated with the second benzene ring that is coupled to the right lead. The torsional 
angle ip of the two benzene rings is explicitly taken into account as a dynamical 
variable, since it modulates the overlap of the 7r-orbitals of the benzene rings and 
this will have an important effect on the conductivity of the molecular junction. 
The influence of the other vibrational modes is neglected in the present study. The 
Hamiltonian of the model thus reads 



The first line of Eq. describes the torsional dynamics of the two benzene rings in 
the absence of an excess electron, the second line describes the electronic resonance 
states corresponding to the situation when the excess electron is located on the 
molecular bridge, and the last line represents the free electrons in the leads and the 
coupling between the leads and the two n— orbitals. The kinetic energy operator 
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for the torsional motion is given by 

where I = ^Md 2 is the moment of inertia, M denotes the mass of the carbon 
atom and d is the diameter of the benzene ring. The potential energy Vo(ip) for the 
torsional motion of the isolated biphenyl molecule can be calculated by standard 
methods of quantum chemistry. In this study, we approximate the potential in form 
of a three-term Fourier series 

V (<p)=C + f[l-cosM], (3) 

i=2,4,6 

where the parameters Ci have been adopted from Ref. (2H| • We furthermore consider 
a symmetric molecular junction with e\ = 62 = eo, Vkeh = ^fcek = ^ k ano - ^fceL = 
^fceR = 0- The electronic coupling between the two rings is assumed to be of the 
form 

p(<p)=pocoB{<p). (4) 

The parameters eo and (3q can in principle be determined from electronic- 
structure calculations of the potential-energy surface of the molecular anion. The 
vertical and adiabatic electron affinities of biphenyl have been calculated, e.g., in 
Ref. [301 yielding slightly different result for different levels of electron structure 
theory. Since the molecular electronic structure is modified by the coupling of the 
biphenyl to the leads, we cannot directly adopt these data, but they can guide us 
to select the parameters in a reasonable range of values. Within our model, the 
adiabatic potential-energy curves of the molecular anion can be obtained by diag- 
onalization of the first two lines of the Hamiltonian 0} for fixed nuclear geometry 
(i.e. T v «0). This way, we obtain two potential energy curves 

Vt\<fi) = V (<p)+eo ± IAjcosMI, (5) 



corresponding to the ground and excited electronic state of the molecular anion, 
respectively. To determine the parameters eo and fto, we impose the following re- 
quirements as motivated by the results of the electronic structure calculations [301 
mentioned above: We require the equilibrium geometry of the anion to be planar, 
i.e. (p = 0, with a small, but positive, adiabatic electron affinity and, furthermore, 
assume a negative vertical electron affinity of the order of -O.leV at the equilibrium 
geometry of the neutral molecule (ipa — 42°). These requirements are fulfilled for 
eo = 0.9eV and (3q = leV. The corresponding potential-energy functions Vo((p) 
(full line) and \ip) (dashed line) are shown in Fig. [21 The upper of the two 
adiabatic molecular-anion potentials, (cp), which exhibits an avoided crossing 
with (<p), is too high in energy to be seen in Fig. [3 
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Fig. 2. Potential-energy curves for the model (for details see text). The energy levels of 
the hindered torsional motion of the neutral molecule are also shown as horizontal lines. 
Note that the states with energies well below and well above zero are nearly degenerate 
(four-fold degenerate below and two- fold degenerate above zero energy, respectively). 



As is well known (see for example Ref. [2]), the effect of the electronic coupling 
to the leads is fully described by specifying the self-energy functions 

Ey(e) = £(&|ff|&>^— (tfkTOi) = <y«Si(e). (6) 
k e ek 

The presence of the factor 5ij is a consequence of the fact that we neglect the direct 
coupling of the right benzene ring to the left lead and the coupling of the left ring 
to the right lead. As in our previous work , we assume that the self-energies are 
proportional to the Hubbard Green's function 



2a 2 



which is the exact self-energy for a one dimensional semi-infinite atomic chain with 
nearest-neighbor tight-binding interaction. Within this model, the energy of an 
electron in state \4>k) in lead j reads 

e k = Hj +20cosk. (8) 

The nearest-neighbor coupling parameter /3 determines the width of the conduction 
band, which is 4/3, and a is the coupling strength between the last atomic site of the 
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lead and the 7r-orbital \<pj) of the adjacent benzene ring. For simplicity we assume 
this coupling to be the same for both leads and independent of the angle ip. The 
chemical potential is different for both leads, the difference eV = /Xl — Mr being 
the bias voltage across the junction. By analytic continuation we obtain the real 
and imaginary parts of the self-energy Ej = Aj — ^Tj on the real energy axis. The 
imaginary part (which is sometimes also called 'energy-dependent decay width of 
the resonance state') is given by 

r ^) = ^\/ 4/32 -( e -^ 2 W 

inside the conduction band (|e — fij\ < 2(3) and Tj = otherwise. In our model 
studies we use a value (3 — leV, which corresponds to a conduction-band width of 
4 eV. The strength a of the coupling between the molecule and the leads is varied 
as described below. 



2.2 Calculation of transmission probability and current 

Employing scattering theory, the transmission probability for scattering of a con- 
duction electron from the initial state \4>ki) in the left lead to the final state \4>k t ) in 
the right lead, accompanied by a vibrational transition of the molecule from initial 
state \Xvi) to final state \xv t ) is given by 

t L ^n(euVue u V{) = 6(E Vf +e i -E v -e i )r 1 {e i )T 2 (e i ) \(xv { \[E -H d - F{E)]^\ Xvi )\ 2 ■ 

(10) 

Thereby, the Dirac delta function accounts for energy conservation 

E = E Vi + ei = E vt +€{, (11) 

with £j (ef) and E Vi (E V{ ) being the initial (final) electron and vibrational energies, 
respectively. H d is a 2 x 2 matrix in the electronic basis {|</>i), |<fe)} given by 

n _ ( T„ + V ftO+eo /%) \ h9) 

d \ m T v + V (<p)+e )> (LZ) 

where each matrix element is a vibrational operator. Similarly, F(E) is a diagonal 
2x2 matrix with matrix elements 

F : (E) ee Ej(E -T v - Vo(v)) = E \Xv)^(E - E v )( Xv \. (13) 

V 

Since the operator T v + Vo(tp) does not commute with 0(<p), the transmission prob- 
ability cannot be obtained in a closed analytical formula. However, the expression 
l|l(J|) can rather easily be evaluated numerically adopting, for example, the free-rotor 
basis for the theoretical treatment of the torsional motion of the bridge. 
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Based on the transmission probability (|1C)|> . the current through the molecular 
bridge is calculated employing the generalized Landauer formula ^3 

I = -^Pvi I de; / de{{f L (cj)[l - f R (ef)]t h ^ R (ei,Vi,ef,Vf) 

Vi,Vf 

-/fl(ei)[l - fLicf^tR^ieuVueuVf)} . (14) 

Here, 

p Vi = exp(-E Vi /kT)/Z (15) 

denotes the Maxwell-Boltzmann distribution of the initial vibrational states for a 
given temperature T and 

fj{E) = i +exp [(E- H )m (16) 

is the Fermi-Dirac distribution for the conduction-electrons in the leads. The bias 
voltage enters the formula through the chemical potentials = eV/2 and /ir = 
—eV/2 for the left and right leads, respectively, and is thus also present in the 
transmission function through the dependence of Tj (e) on the chemical potentials 
Mr 

3 Results and discussion 
3.1 Transmission probability 

The transmission probabilities tL->ii(ei, i>i,e{,vt) for zero bias have been calculated 
according to formula IjlOl) employing a free rotor basis with typically 200 basis 
functions. To illustrate the results, we average £l— >R over the lowest four vibrational 
energy levels V{ — 0, 1, 2, 3 (which are nearly degenerate due to the symmetry of the 
torsional potential), sum over all final states Vf and integrate over the final energy 
of the electron ef. The resulting transmission probability is shown for three values 
of the coupling strength (a = 0.05 eV, 0.1 eV and 0.2 eV) in figure It is seen 
that the transmission is essentially nonzero only in the region where the molecular 
resonance is located. The process can thus be understood as resonant tunneling of 
the electron from the left lead to the right lead through the eigenstates of These 
states coincide (in the relevant energy region) with the torsional eigenstates in the 
potential V} \tp) (shown as dashed line in Fig.[3J). In the transmission probability, 
the individual states are resolved only for the smallest value of the coupling a = 0.05 
eV. For larger values of a, the individual peaks are wiped out due to the stronger 
coupling with the leads. 

It is also interesting to compare the full transmission probability with the elastic 
one, calculated for frozen vibrations 

&V(e u e uV ) = ( 5(e f -e i )r 1 (e i )r 2 ( ef ) y . 

(17) 
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Fig. 3. Energy dependence of the averaged transmission probability (see text) for three 
coupling strengths a and temperature T = 0. The adiabatic nuclei approximation (dashed 
lines) is compared with the full nonadiabatic calculation (full lines). The elastic transmis- 
sion calculated at the equilibrium geometry of the neutral molecule is also shown for one 

case (dotted line). 

The elastic transmission function depends on the instantaneous value of the 

torsional angle ip. Fig. [21 shows this function (dotted line) at the equilibrium ge- 
ometry of the neutral molecule (ipo = 42°), integrated over the final energy of the 
electron e/ for one value of the electronic coupling (a = 0.1 eV). It is seen that 
the elastic transmission is much narrower than the full inelastic transmission, thus 
demonstrating the pronounced broadening effect of the torsional degree of freedom. 
To investigate if this is a truly dynamical effect or rather the result of the distribu- 
tion of torsional angles which contribute to the quantum-mechanical initial state 
IXuj), we have averaged £l-2-r ovcr the initial distribution of the torsional angle ip 
for temperature T. Integrating furthermore over the final energy of the electron, 
we obtain the expression 



which in the context of electron-molecule scattering is also called the adiabatic- 
nuclei approximation . As the results in Fig.|21show, the adiabatic- nuclei approx- 
imation works very well for the system under consideration except for extremely 
small coupling to the leads where individual vibrational states of the molecular 
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anion are resolved, which is not reproduced in adiabatic-nuclei approximation. 
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Fig. 4. Dependence of the transmission probability on the final torsional state for three 
different initial energies of the electron: below resonance (ti = 0), on resonance (ej = 0.15 
eV) and above resonance (a = 0.4 eV). The vertical dotted line shows the position of the 
state with energy E v — (onset of freely rotating states) (cf. Fig. yj. 

Being based solely on elastic transmission mechanisms, however, the adiabatic- 
nuclei approximation cannot describe the vibrational excitation of the molecular 
bridge accompanying the transmission of the electron. This information, which is 
important to characterize the possible heating of the molecular bridge, is contained 
in the full transmission probability (|10|l . As an illustration, Fig.HJshows the depen- 
dence of the inelastic transmission probability (integrated over the final electron 
energy) on the final torsional state for a coupling strength of a = 0.1 eV and three 
different initial energies of the electron. The results exhibit a substantial vibrational 
excitation of the bridge in the case of resonant electron transmission. The effect 
of the electronic coupling strength a on this resonant heating process is studied in 
Fig. It is seen that highly excited torsional states are not populated for strong 
coupling to the leads (corresponding to a broad resonance), although in this case 
the total transmission probability is larger than for small values of a (cf . Fig. EJ) . 
This result is due to the fact that for strong coupling between molecular bridge and 
leads the residence time of the electron on the molecular bridge is too short to allow 
efficient electron-nuclear interaction and thus the total transmission is dominated 
by the elastic contribution. 

In the context of electron transport through molecular bridges, vibrational ex- 
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Fig. 5. Dependence of the transmission probability on the final torsional state for three 
different strengths of the coupling of the molecule to the leads and an initial energy of the 

electron, ei = 0.15 eV. 



citation (and possibly dissociation) of the molecular bridge is presumably the most 
important process induced by the coupling between electron and vibrational degrees 
of freedom. Another interesting process is electron assisted tunneling between the 
wells of the torsional potential Vo(ip). To study this process, we consider torsional 
states which are localized in the wells of the potential Vo(ip). The potential Vo((p) 
has four equivalent equilibrium geometries (po,i = ±</?o an d ^2,3 = ^ ± Vo (see Fig. 
[2J). As a consequence, combinations of the four nearly degenerate states with lowest 
energy can be used to obtain states which are localized in the four wells. Due to 
the near-degeneracy, these states are close to stationary states. In Fig. [5] we show 
the transmission probability of the electron starting from a torsional state located 
at tpo and ending up in the same state (upper two curves) or in a torsional state 
that is localized in the adjacent potential well at (pi (lower two curves). It is seen 
that the resonance peak for both processes is located at the same position, the 
width of the peak, however, is significantly narrower for the lower two curves, i. e. 
for current-induced tunneling between the potential wells. This can be understood 
as suppression of the electron transmission due to (assisted) tunneling through the 
torsional barrier. Another interesting result is that the transmission probabilities 
for the two different processes at energies close to the resonance peak become very 
similar for smaller coupling (a = 0.05 eV, full lines). This is a consequence of the 
long-lived character of the molecular anion in the case of small coupling to the 
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Fig. 6. Electron-assisted tunneling between adjacent wells in the torsional potential. 
Shown is the energy dependence of the electron transmission probability starting from a 
torsional state which is localized around tpo with (lower two curves) and without (upper 
two curves) transition to another, equivalent torsional state in the neighboring well. 



leads. For sufficiently long lifetimes of the molecular anion, the torsional state in 
the potential vi ' has the same probability to decay to both neighboring wells in 
the neutral potential Vq. 

3.2 Current-voltage characteristic 

Finally, we consider the current-voltage characteristic for our model. To obtain the 
current for a given voltage via Eq. (|14|l . we have calculated the transmission coeffi- 
cients (f 101) for a series of bias voltages on a dense energy grid. Fig. [7] shows results 
for different temperatures T. The I — V characteristics exhibit a step-like structure. 
Due to thermal broadening, the width of the step increases with increasing tem- 
perature, ft should be emphasized, that this thermal broadening is mainly a result 
of the thermal excitation of the torsional degree of freedom (i.e. Eq. I|15|l ). whereas 
the thermal effect of the electrons (in Eq. I|16[l ) is almost negligible. 

The adiabatic-nuclei approximation (data not shown) follows the full calcula- 
tion closely. For voltages above 0.3 V, the difference is not noticeable. It becomes 
more apparent close to zero bias. For example, for T ~ 100 K the zero-bias conduc- 
tivity is 4.5 x 10 -4 , whereas the adiabatic nuclei approximation gives 6.7 x 10~ 4 . 
This difference is a consequence of the energy loss during the transmission process 
which is correctly described in the full calculation, but not in the adiabatic-nuclei 
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Fig. 7. Current-voltage characteristics for different temperatures (see legend) and an 
electronic coupling strength a = O.leV. 



approximation. 



4 Conclusions 



We have analyzed vibrationally inelastic electron transport through a molec- 
ular junction in the presence of an anharmonic large-amplitude motion of the 
molecular bridge. The study was based on a mechanistic model for a biphenyl 
molecule between two metal electrodes. The results of the study demonstrate the 
effects of electron- vibrational coupling on the transmission probability and the cur- 
rent through the junction. While for small coupling between the molecule and 
the leads the torsional degree of freedom gives rise to well resolved structures in 
the transmission function, for stronger coupling it mainly results in a broadening, 
which increases significantly for higher temperature. We have, furthermore, studied 
the vibrational excitation process accompanying the transmission of the electron 
through the molecular bridge, which is particularly pronounced for electron ener- 
gies close to the molecular resonance state and/or small to moderate electronic 
coupling strength. A correct description of the vibrational excitation process is also 
important to characterize the heating of the molecule caused by the electron trans- 
port. To investigate this heating process in more detail, the coupling to the other 
vibrational degrees of freedom of the molecule (which were neglected in the present 
study) and the environment has to be taken into account [211 1221 123) . 
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Another interesting process, which is not related directly to vibrational exci- 
tation of the molecular bridge, are structural changes of the bridge, in particular 
electron assisted tunneling between different wells of the torsional potential. Our 
studies show that electron assisted tunneling becomes particularly important when 
narrow resonances are present. 

We have also investigated the validity of approximate methods, in particular 
the purely elastic description of the transmission process and the adiabatic nuclei 
approximation. As was found in our previous study for models with small ampli- 
tude motion |23| . the purely elastic calculation typically predicts much too narrow 
transmission probabilities. The adiabatic nuclear approximation, on the other hand, 
which takes into account the 'static' effect of the vibrational degrees of freedom due 
to their initial distribution, can describe the average transmission probability rather 
well, but it is not capable of describing the vibrational excitation accompanying the 
electron transport. Furthermore, it is also expected to fail if the density of states 
of the leads varies quickly in the energy range of interest. This aspect, which is 
expected to be of particular importance if the molecular bridge is bound to semi- 
conductor instead of metal electrodes will be the subject of future work. 
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